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Abstract 

^ Strong interactions between hadrons include attractive and repulsive components. In statistical 

approaches attractive interactions are modelled by including resonances while repulsive interactions 
by taking into account a finite proper volume of particles. In this paper the influence of repulsive 
interactions on matter properties is considered within the excluded volume van der Waals hadron- 

a 

resonance gas. Quantitative results are presented for matter at the chemical freeze-out in central 
nucleus-nucleus collisions at relativistic energies. In particular, it is shown that repulsive interactions 
lead to a significant decrease of collision energy at which the net-baryon density has a maximum. 
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I. INTRODUCTION 



Statistical models of the hadron gas are an important tool to extract properties of matter 
created in relativistic nucleus- nucleus collisions (see, e.g., Refs. [IHZ1)- Basic parameters of 
these models are the matter temperature T, baryon chemical potential and volume V. If 
supplemented by additional model parameters monitoring deviations from the chemical equi- 
librium [8], they approximately fit rich data on mean hadron multiplicities in a broad range of 
reactions, from e + +e~, p+p, and p+p [9] at low energies to central Pb+Pb collisions at the 
highest LHC energy [ID] . 

The most popular version of the statistical models of hadron matter is the ideal hadron- 
resonance gas (I-HRG), i.e., a statistical system of non-interacting hadrons and resonances. 
One argues, based on the Dashen, Ma and Bernstein theorem [IT], that resonances introduced 
to the ideal hadron gas take into account attractive interactions between hadrons. The repulsive 
part of interactions between hadrons is usually accounted for by the van der Waals excluded 
volume procedure generalized to the relativistic case of a variable number of hadrons [12] . The 
resulting excluded volume model is not anymore an ideal gas model, and in this paper it will 
be denoted as the EV-HRG model. Another popular example of modelling attractive and 
repulsive interactions between hadrons is the relativistic mean field theory in a form of the 
Walecka model [13] and its different modifications (see e.g., recent paper [14J and references 
therein). In this approach, scalar and vector meson fields describe respectively the attractive 
and repulsive forces between baryons. 

Both attractive and repulsive interactions are important for qualitative as well as quantita- 
tive description of properties of strongly interacting matter. For example, the nucleon-nucleon 
potential includes both parts - attractive at large and repulsive at small distances. The pres- 
ence of both attractive and repulsive interactions between nucleons is crucial for the existence 
of stable nuclei. Moreover an important undesired feature of the I-HRG model at high tem- 
peratures was noted in Refs. [15j [IB] . Due to the large number of different types of baryons 
and mesons, the point-like hadrons would always become the dominant phase at very high 
energy density. Just the excluded volume effects ensure a transition from a gas of hadrons and 
resonances to the quark-gluon plasma. 
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The aim of this paper is to recall the role of repulsive interactions between hadrons. It 
is shown that the excluded volume hadron-resonance gas model yields different properties of 
matter at the freeze-out than the ideal hadron-resonance gas, if densities and their collision 
energy dependence are considered. Quantitative results are presented for two examples, namely 
collision energy dependence of the net-baryon density [T7] as well as the ratio of baryon and 
meson entropy density [TSH2T)] . They are selected, because of conjectures that the maximum 
baryon density as well as the transition between baryon and meson dominated matter may 
be related to the onset of deconfmement observed in central Pb+Pb collisions at the CERN 
SPS energies [21] . Note that chemical freeze-out parameters T and //^ in nucleus-nucleus 
collisions are straightforwardly connected to the data on hadron multiplicities. The energy 
range considered in this paper is presently studied experimentally at the CERN SPS [22J and 
the BNL RHIC [23]. In future this effort will be extended by experiments at new accelerators, 
at the JINR NICA [24] and the FAIR SIS-100 [251 [26]. 

The paper is organized as follows. In Sec. [TT] the ideal hadron-resonance gas model is 



introduced. The excluded volume hadron-resonance gas model is presented in Sec. Ill, where 
also quantitative results for densities at the chemical freeze-out in central Pb+Pb collisions are 



given and discussed. A summary given in Sec. IV closes the paper. 



II. IDEAL HADRON RESONANCE GAS 



In the grand canonical ensemble the pressure of the I-HRG is given by 
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(1) 



where T is the system temperature, 77 = — 1 and rj = 1 for bosons and fermions, respectively, 
while rj = corresponds to the Boltzmann approximation. For a hadron i, rrii is its mass and 
di is the spin degeneracy. The chemical potential is given by 



Hi = hjJB + SiHs + QiHQ 



(2) 
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with hi = 0, ±1, Si = 0, =Fl, =f2, =)=3 and qi = 0, ±1, ±2 for hadrons. The number density of a 
hadron i reads: 

-l 



<9/lj 27T 



, -Hi , , 

exp | | ± 1 



(3) 



Considering the temperature T, baryon chemical potential \Lb, and volume V as free parameters 
one can fit mean hadron multiplicities (Ni) = Vrii measured in relativistic nucleus-nucleus 
collisions at each collision energy. In this analysis, /is and //q are expressed as functions 
of T and \xb when the conditions on strangeness, (S) = 0, and electric to baryon charge 
ratio, (Q)/{B) = Z/A, are taken into account. Most of experimental data on nucleus-nucleus 
collisions concern yields of long-lived hadrons, which include products of resonance decays. 
This requires a proper treatment of short lived resonances, namely products of their strong and 
electromagnetic decays should be added to mean multiplicities of stable hadrons. In this paper 
the numerical implementation of the hadron-resonance gas model provided by the THERMUS 
package [27] is used to calculate the relevant quantities according to Eqs. (|4]{6]). Particles and 
resonances (all mesons up to _fQ(2045)) and baryons (up to quantum statistics as well as 
width of resonances are included. 

The analysis of central Pb+Pb (Au+Au) collisions registered by experiments at SIS, AGS, 
SPS, and RHIC allows to establish the collision energy dependence of T and fiB which can be 
parameterized as [T9] : 

T = 166 - 139 u% - 53 u% , u B = , (4) 

^ B ^ B ' ^ B 1000 + 0.273 t/s]^ ' V ; 

where y/s^jy is the center of mass energy of a nucleon pair. Values of T, and a/s/vjv in 
Eq. Q should be given in MeV. The chemical freeze-out line, T = T(/i#), as well as the energy 
dependence of the T and /j,b parameters are shown in Fig. [l] (a) and (b), respectively. 

The net-baryon density ps, entropy density s and energy density e can be found from the 
system pressure p using the thermodynamical relations: 

dp dp dp dp 

PB = Wb ' S = df ' £ = T df + ^- p - (5) 

With chemical freeze-out parameters given by Eq. (|4]) and the ideal gas expression Eq. Q 
for the system pressure, one finds the quantities in Eq. ^ as functions of collision energy. 




The I-HRG model is based on the assumption of complete thermal and chemical equilibrium. 
An additional I-HRG parameter, the strangeness suppression factor 75, has to be introduced 
to account for deviations of strange hadron multiplicities from chemical equilibrium [8]. Its 
dependence on T and fiB obtained by fitting hadron yields measured in the full phase-space 
can be parameterized as [7] 

75 = 1 -0.396 exp (- h23 jP) • ( 6 ) 

At the AGS and SPS energies the 75 parameter is significantly smaller than one, which means 
under- saturation of strange hadron yields with respect to the chemical equilibrium. The 75 
parameter should be included in the model if a proper description of strange hadron yields 
is required. However, the relation between T and /ig as well as the chemical freeze-out line 
Eq. Q, obtained within the models with and without the 75 parameter is similar. 

The baryon density p 1 ^ as a function of collision energy calculated following the freeze-out 
line Q is shown in Fig. [2] (a) . In this and the following figures the laboratory collision energy 
per projectile nucleon E iah is used to present the dependence on collision energy. Its connection 
to center-of-mass energy per nucleon pair, ^/snn, is given by ^Snn = ^J^rriN Ei a b + 2m 2 N , 

5 



0.14 - 



PQ 



0.07 - 




E lab [AGeV] 



H 

C/3 9 



6 - 



\ 

Baryons + Mesons 


4=1 

g, < 1 , Boltz 


— ^-^Sl 






(b) 



10 100 

E, ab [AGeV] 



FIG. 2: The net-baryon density ps (a) and the ratio s/T 3 (b) along the chemical freeze-out line 
Eq. Q and 75 according to Eq. ^ are shown by the solid lines. The dashed lines correspond to 
75 = 1 for pb (a) and s/T 3 (b). The dotted line corresponds to the Boltzmann approximation in 
s/T 3 . 



where is the nucleon mass. As seen in Fig. [2] (a), the baryon density has a maximum [T7j 
at Eiaf, = 34A GeV. This is the collision energy at which the NA49 Collaboration observed the 
maximum of the K + /n + ratio (the horn) and other signals of the onset of deconfmement [21] . 

The total entropy density as a function of collision energy following the freeze-out line 
Eq. Q is shown in Fig. [2] (b). Meson sm and baryon sb entropy densities are also presented in 
the figure. With increasing collision energy the baryon-dominated (sb > sm) matter changes 
to meson-dominated (% > sb) matter. For the I-HRG model this transition is located at 
E lab = AQA GeV. 

For the T—[Lb values at the chemical freeze-out line Eq. Q the role of quantum statistics 
is small. For baryons the Fermi statistics changes their densities by less than 1%. The largest 
density change due to the Bose statistics is for pions. It is, however, still smaller than 10%. The 
ratio s/T 3 calculated within the Boltzmann approximation, i.e. 77 = in Eq. ([!]), is shown in 
Fig. [2] (b) by the dotted line. The deviations from the results with quantum statistics included 
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are hardly visible. They are even smaller for ps and thus the corresponding line calculated 
with the quantum statistics included is not plotted. The collision energy dependence of ps and 
s/T 3 calculated for the 75 parametrization Eq. ^ and for 75 = 1 is also shown in Fig. [2J One 
concludes that the energy at which p# has the maximum as well as the energy of the transition 
between baryon-dominated and meson-dominated matter are approximately independent of 
the quantum statistics and the degree of strangeness equilibration. 

It was suggested [H] that the maximum of net-baryon density and/or the transition from 
baryon to meson dominance may be related to the anomalous behavior of the K + /tc + ratio |21j. 
In the next section these phenomena are examined by taking into account repulsive interactions 
between hadrons. 



III. EXCLUDED VOLUME HADRON RESONANCE GAS 

The results presented in Section II have been obtained within the ideal hadron-resonance 
gas model in which only attractive interactions between hadrons are taken into account by the 
inclusion of resonances. In this section the role of repulsive interactions is considered within 
the excluded volume hadron-resonance gas model. 

The van der Waals excluded volume procedure corresponds to a substitution of the system 
volume V by the available volume V — Ei^i-^i , where Vi = 4(4/37rrf) is the excluded volume 
parameter and rj is the corresponding hard sphere radius of a particle i. This result, in partic- 
ular, the presence of a factor of 4 in the expression for Vi, can be rigorously obtained for a low 
density gas of particles of a single type (see, e.g., Ref. [28]). In the grand canonical ensemble, 
the substitution V — > V — Ej t^A^ leads to a transcendental equation for the pressure of the 
EV-HRG 1 [HI US]: 

P = ^2pf( T ,fii) ; fii = Pi - ViP . (7) 

i 

Using Eq. Q one finds the net-baryon, entropy and energy densities: 

_ EjMgg^j s T,i _ T^£iM ( 8 ) 

PB 1 + Ej vjnf(T, fit) ' 1 + Ei v j n f(T, fij) ' 1 + Ei vjn?(T, fij) ' 1 ) 

1 A discussion of other excluded volume formulations can be found in Ref. 1291 . 
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In comparison to the corresponding densities calculated within the I-HRG model the densities 
in the EV-HRG model Eq. ^ are lower because of two reasons: 

i) due to the suppression factor [1 + £\ Vjii^fT, /I,)] -1 and 

ii) due to the shift in chemical potential (ii — » fii which in the Boltzmann approximation 
leads to the suppression factor exp(— Vip/T) < 1. 

The shift of the chemical potential makes the Boltzmann approximation even more accurate 
than in the case of the ideal gas. If all proper volume parameters are the same v i = v (i.e. 
V{ = r), the Boltzmann approximation gives the total suppression factor R 

R{T,fi B ; r) = — , 9 

1 + v J2j n) d {T, Hj) 

the same for all densities of Eq. ([8]): 

p B (T,fi B ) = Rp^(T,p B ) , s(T,fi B ) = Rs id (T,p B ) , e(T,p B ) = Rp%(T,fi B ) (10) 

and rii(T, p B ) = Rnf(T, p B ) . Typical values of hard-core radii considered in the literature j2j 
130] are r = (0.3 -r- 0.8) fm. The energy dependence of the suppression factor Eq. ^ calculated 
along the chemical freeze-out line for r = 0.5 fm and r = 1 fm is shown in Fig. [3} The R factor 
(|9| decreases monotonously with increasing collision energy. For example, for r = 0.5 fm one 
finds R = 0.9 and R = 0.4 at small and large Ei^, respectively. One may therefore expect a 
decrease of the value of p B at its maximum by a factor of 0.5 (for r = 0.5 fm), and a shift of 
the position of the maximum to a smaller collision energy. In fact, in Fig. [4] (a) one observes 
that the maximum of the net-baryon density is located at Ei a b = 17 A GeV for r = 0.5 fm 
and at E lab = 7A GeV for r = 1 fm, instead of E lab = 3AA GeV for the I-HRG model. It is 
also seen that the value of p B at the maximum decreases strongly with increasing value of the 
hard-core radius. The entropy density shown in Fig. [4] (b) is reduced by the same suppression 
factor. The collision energy at which baryon and meson entropy densities are equal is, however, 
independent of R and is located at Ei a b = 46A GeV. This is true only if the hard-core radius 
parameter is the same for all hadrons. 

It is instructive to consider the role of the excluded volume effects for different hard-core 
radii of baryons r B and mesons tm- As an example, the results for r B = 0.5 fm and tm = are 
presented in Fig. |4| In a comparison to the results for r = 0.5 fm for all hadrons one observes 
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FIG. 3: The excluded volume suppression factor R Eq. Q as a function of Ei ab along the chemical 

The solid and dotted lines correspond to r = 0.5 fm and r = 1 fm, 
respectively. 

small changes of but a significant shift of the transition point between baryon and meson 
dominated matter. Its position decreases from Ei ab = 46A GeV to E [ab = 23A GeV. The model 
with non-equal radii of hadrons needs, however, further detailed studies. This is because fits 
to the hadron yields performed with the EV-HRG model with non-equal radii should give 
different freeze-out parameters T and [is than those in Eq. Q obtained within the I-HRG 
model. This analysis is beyond the scope of the paper. 

IV. SUMMARY 

The ideal hadron-resonance gas model in the grand canonical ensemble is simple and has 
only few free parameters. In spite of this it is successful in describing the bulk properties of 
mean hadron multiplicities in high energy collisions. The model takes into account attractive 
interactions between hadrons but ignores repulsive interactions. The latter are needed to catch 
basic qualitative features of strong interactions, e.g, the phase transition between hadron- 
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FIG. 4: (a): The baryon density along the chemical freeze-out line from Eqs. (4j6). Dashed-dotted 
line corresponds to the model with r = 0, dashed line to = 0.5 fm and r^j = 0, solid line to 
r B = r M = r = 0.5 fm, and dotted line to rg = tm = r = 1 fm. (b) The ratios sb/T 3 and sm/T 3 
along the chemical freeze-out line Eqs. (4j6). Dashed-dotted lines corresponds to the model with 
r = 0, solid lines to r& = tm = r = 0.5 fm, and dashed-dotted lines to rs = 0.5 fm and tm = 0. 



resonance gas and the quark-gluon plasma. Moreover, they strongly modify properties of the 
hadron-resonance gas. The most common way to include repulsive interactions in the hadron- 
resonance gas model is to follow the van der Waals excluded volume procedure and introduce 
hard-core radii of hadrons. 

If radii of all hadrons are assumed to be the same, the chemical freeze-out parameters, 
temperature and baryon chemical potential, fitted to data on mean hadron multiplicities are 
identical to those obtained within the ideal hadron-resonance gas model. However, all densities 
calculated within the van der Waals model are lower than the corresponding densities obtained 
within the ideal gas model and thus the fitted volume parameter in the van der Waals gas 
formulation is significantly larger. The density suppression factor R depends on the T and 
[is parameters, which in turn depend on collision energy. Consequently the collision energy 
dependence of densities is sensitive to the assumed hard-core radius of hadrons. In particular, 
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the energy at which net-baryon density has a maximum decreases from about Ei ab = 34A GeV 
for the ideal gas model to about Ei ab = 7 A GeV for the excluded volume model with r = 1 fm. 

If the radii of hadrons are assumed to be different, densities of different hadrons are modified 
differently. In particular, the collision energy at which the transition between baryon and 
meson dominated matter takes place decreases from Ei ab = 46A GeV for the ideal gas model 
to E lab = 23A GeV for the excluded volume model with the hard-core radii of baryons and 
mesons equal to 0.5 fm and 0, respectively. Clearly, the excluded volume effects are even larger 
for the hadron matter at stages preceding the chemical freeze-out in nucleus-nucleus collisions, 
i.e. at larger values of the energy density. 

In view of this study, the suggestions that the collision energy of the onset of deconfinement 
coincides with the energy at which the net-baryon density at the chemical freeze-out reaches its 
maximum and/or the transition between baryon and meson dominated matter takes place are 
premature. This is because the results obtained within the hadron gas approach are strongly 
dependent on the modelling of repulsive interactions between hadrons. Experimental and/or 
theoretical methods to better estimate hard-core radii of hadrons within the excluded volume 
model are needed to improve our understanding of the properties of hadron-resonance matter. 
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